Quantum-to-classical crossover of mesoscopic conductance fluctuations 
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We calculate the system-size-over-wave-length (M) dependence of sample-to-sample conductance 
fluctuations, using the open kicked rotator to model chaotic scattering in a ballistic quantum dot 
coupled by two TV-mode point contacts to electron reservoirs. Both a fully quantum mechanical 
and a semiclassical calculation are presented, and found to be in good agreement. The mean 
squared conductance fluctuations reach the universal quantum limit of random-matrix-theory for 
small systems. For large systems they increase oc M 2 at fixed mean dwell time td oc M/N. The 
universal quantum fluctuations dominate over the nonuniversal classical fluctuations if TV < y/~M. 
When expressed as a ratio of time scales, the quantum-to-classical crossover is governed by the ratio 
of Ehrenfest time and ergodic time. 

PACS numbers: 73.23.-b, 73.63.Kv, 05.45.Mt, 05.45.Pq 



I. INTRODUCTION 

Sample-to-sample fluctuations of the conductance of 
disordered systems have a universal regime, in which they 
are independent of the mean conductance. The require- 
ment for these universal conductance fluctuations 1 * 2 - is 
that the sample size should be small compared to the lo- 
calization length. The mean conductance is then much 
larger than the conductance quantum e fh. 

The same condition applies to the universality of 
conductance fluctuations in ballistic chaotic quantum 
dots^, although there is no localization in these sys- 
tems. Random-matrix-theory (RMT) has the universal 
limit 



lim var G 

N^oo 



1 



(1.1) 



for the variance of the conductance G in units of e 2 fh. 
Here iV is the number of modes transmitted through each 
of the two ballistic point contacts that connect the quan- 
tum dot to electron reservoirs. Since the mean conduc- 
tance (G) = N/2, the condition for universality remains 
that the mean conductance should be large compared to 
the conductance quantum. 

In the present paper we will show that there is actually 
an upper limit on TV, beyond which RMT breaks down 
in a quantum dot and the universality of the conduc- 
tance fluctuations is lost. Since the width W of a point 
contact should be small compared to the size L of the 
quantum dot, in order to have chaotic scattering, a triv- 
ial requirement is N -C M, where M is the number of 
transverse modes in a cross-section of the quantum dot. 
(In two dimensions, N ~ W/Xp and M ~ L/Xp, with 
Af the Fermi wavelength.) By considering the quantum- 
to-classical crossover, we arrive at the more stringent re- 
quirement 



l«f« \/Me ATc ' s/2 , 



(1.2) 



more stringent than N <§C M because, typically, A -1 and 
T crg are both equal to the time of flight to across the 
system, so the exponential factor in Eq. (|1.2|) is not far 
from unity. 

Expressed in terms of time scales, the upper limit in 
Eq. I|1.2I) says that r crg should be larger than the Ehren- 
fest time 5-6 



te = max 



0,A -1 ]ii 



N 2 
~M 



(1.3) 



with A the Lyapunov exponent and T erg the ergodic time 
of the classical chaotic dynamics. The requirement is 



The condition r org > te which we find for the universality 
of conductance fluctuations is much more stringent than 
the condition td > t~e for the validity of RMT found in 
other contexts.SiWiMiWi 11 ^ 2 ! 1 ? Here t d « (M/N)t is 
the mean dwell time in the quantum dot, which is ^> T crg 
in any chaotic system. 

The outline of this paper is as follows. In Sec. II we de- 
scribe the quantum mechanical model that we use to cal- 
culate var G numerically, which is the same stroboscopic 
model used in previous investigations of the Ehrenfest 
tim oW 4 . The data is interpreted semiclassically in Sec. 
Ill, leading to the crossover criterion (|1.2[) . We conclude 
in Sec. IV. 



II. STROBOSCOPIC MODEL 

The physical system we have in mind is a ballistic 
(clean) quantum dot in a two-dimensional electron gas, 
connected by two ballistic leads to electron reservoirs. 
While the phase space of this system is four-dimensional, 
it can be reduced to two dimensions on a Poincare surface 
of sectioni£J£ The open kicked rotator 9 ! 11 ! 14 . 17 ! 18 ! 19 ! 20 
is a stroboscopic model with a two-dimensional phase 
space. We summarize how this model is constructed, fol- 
lowing Ref. 0. 

One starts from the closed system (without the leads) . 
The kicked rotator describes a particle moving along a 
circle, kicked periodically at time intervals tq. We set 
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FIG. 1: Classical phase space of the open kicked rotator. The 
dashed lines indicate the two leads (shown for the case td = 
5). Inside each lead we plot the initial and final coordinates 
of trajectories which are transmitted from the left to the right 
lead after at most 3 iterations (with K = 7.5). The points 
cluster along narrow "transmission bands". 

to unity the stroboscopic time To and the Plank con- 
stant H. The stroboscopic time evolution of a wave func- 
tion is given by the Floquet operator T ', which can be 
represented by an M x M unitary symmetric matrix. 
The even integer M defines the effective Planck constant 
h c s = 1/M. In the discrete coordinate representation 
(x m — m/M, m = 0, 1, . . . , M — 1) the matrix elements 
of T are given by 

•Fm'm = M-V 3 e- i */*e a " it3 (-*'» , ' Xm \ (2.1) 

where S is the map generating function, 

S(x',x) = \{x' ~xf - (K/8ir 2 )(cos2irx' + cos27ra), 

and K is the kicking strength. 

The eigenvalues exp(— ie m ) of J- define the quasi- 
energies e m € (0,2-71"). The mean spacing 2ir/M of the 
quasi-energies plays the role of the mean level spacing 8 
in the quantum dot. 

To model a pair of N-mode ballistic leads, we im- 
pose open boundary conditions in a subspace of Hilbcrt 
space represented by the indices m[ . The subscript 
n = 1,2, ...N labels the modes and the superscript 
a = 1,2 labels the leads. A 2N x M projection ma- 
trix P describes the coupling to the ballistic leads. Its 
elements are 

p _/lifm = rc.e {m^}, ^ g\ 

nm \ otherwise. 

The mean dwell time is td — M/2N (in units of tq). 




FIG. 2: Variance of the conductance fluctuations obtained 
numerically by varying e with fixed lead positions. Error bars 
indicate the scatter of values obtained for different lead posi- 
tions. Results are shown as a function of l//i e ff = M, for two 
values of the dwell time td = M/2N. The dashed line is the 
RMT prediction varG = |. 

The matrices P and T together determine the quasi- 
energy dependent scattering matrix 

S(e) = P[e~ %£ - T{1 - P T P)]- 1 TP T . (2.4) 

The symmetry of T ensures that S is also symmetric, 
as it should be in the presence of time-reversal symme- 
try. By grouping together the N indices belonging to the 
same lead, the 2N x 2N matrix S can be decomposed 
into 4 sub-blocks containing the N x N transmission and 
reflection matrices, 

s = ( 2 - 5 ) 

The conductance G (in units of e 2 /h) follows from the 
Landauer formula 

G = Trtt t . (2.6) 

The open quantum kicked rotator has a classical limit, 
described by a map on the torus {x,p \ modulo 1}. The 
classical phase space, including the leads, is shown in Fig. 
n The map relates x,p at time k to x' ,p' at time k+ 1: 

p' = ^S(x',x), p=-^-S(x',x). (2.7) 

The classical mechanics becomes fully chaotic for K > 7, 
with Lyapunov exponent A w h\(K/2). For smaller K the 
phase space is mixed, containing both regions of chaotic 
and of regular motion. We will restrict ourselves to the 
fully chaotic regime in this paper. 

III. NUMERICAL RESULTS 

To calculate the conductance l|2.6|l we need to invert 
the M x M matrix between square brackets in Eq. (|2.4|l . 
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We do this numerically using an iterative procedureJA 
The iteration can be done efficiently using the fast- 
Fourier-transform algorithm to calculate the application 
of T to a vector. The time required to calculate S scales 
as M 2 lnM, which for large M is quicker than the M 3 
scaling of a direct inversion. The memory requirements 
scale as M, because we need not store the full scattering 
matrix to obtain the conductance. 

We distinguish two types of mesoscopic fluctuations in 
the conductance. The first type appears upon varying the 
quasi-energy e for a given scattering matrix S(e). Since 
these fluctuations have no classical analogue (the clas- 
sical map (|2.7|) being e-independent), we refer to them 
as quantum fluctuations. The second type appears upon 
varying the position of the leads, so these involve vari- 
ation of the scattering matrix at fixed e. We refer to 
them as sample-to-sample fluctuations. They have both 
a quantum mechanical component and a classical ana- 
logue. One could introduce a third type of fluctuations, 
involving both variation of e and of the lead positions. 
We have found (as expected) that these are statistically 
equivalent to the sample-to-sample fluctuations at fixed 
e, so we need not distinguish between fluctuations of type 
two and three. 
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FIG. 4: Variance of the classical fluctuations of the transmis- 
sion probability Pi^2 upon changes of lead positions, calcu- 
lated numerically from the map 12.71 . The data is shown for 
four values of the dwell time td, as a function of the Lyapunov 
exponent A = \n(K/2). The dotted lines are the analytical 
prediction 14.31 . with fit parameters c — 1.6 and r crg = 0.68 
(the same for all data sets) . 



IV. INTERPRETATION 
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FIG. 3: Same as Fig. H but now for an ensemble in which 
the lead positions and the quasi-energy are both varied. The 
dashed lines are the sum of the RMT value 14.11 and the 
classical result 14.21 . Results are shown for three values of 
the kicking strength K. Open symbols are for the dwell time 
td = 10 and closed ones for td = 20. 

We have calculated the variance varG = (G 2 ) — (G) 2 
of the conductance, either by varying s at fixed lead 
positions (quantum fluctuations) or by varying both e 
and lead positions (sample-to-sample fluctuations). Since 
the quantum interference pattern is completely different 
only for energy variations of order of the Thouless energy 
1/tb, we choose a number td of equally spaced values 
of e in the interval (0, 2ir). We take 10 different lead po- 
sitions, randomly located at the s-axis in Fig. ^ To in- 
vestigate the quantum-to-classical crossover, we change 
/i c ff = 1/M while keeping the dwell time td — M/2N 
constant. The results are plotted in Figs. |21 and 



We interpret the numerical data by assuming that the 
variance of the conductance is the sum of two contri- 
butions: a universal quantum mechanical contribution 
Vrmt given by random-matrix theory and a nonuniver- 
sal quasiclassical contribution V c i determined by sample- 
to-sample fluctuations in the classical transmission prob- 
abilities. 

The RMT contribution equals 3 ' 4 



MT 



(4.1) 



in the presence of time-reversal symmetry. The classical 
contribution is calculated from the classical map 1)2. 7fl . by 
determining the probability P\^2 of a particle injected 
randomly through lead 1 to escape via lead 2. Since the 
conductance is given semiclassically by G c i = NP\^,2, we 
obtain 



V c 



iV 2 varPi^ 2 - 



(4.2) 



We plot var G = Vrmt + V c i in Fig. (dashed curves), 
for comparison with the results of our full quantum me- 
chanical calculation. The agreement is excellent. 

We now would like to investigate what ratio of time 
scales governs the crossover from quantum to classical 
fluctuations. 

To estimate the magnitude of the sample-to-sample 
fluctuations in the classical transmission probability, we 
use results from Ref. |6|. There it was found that the 
starting points (and end points) of transmitted trajecto- 
ries are not homogeneously distributed in phase space. 
Instead, they cluster together in nearly parallel, narrow 
bands. These transmission bands are clearly visible in 
Fig. n The largest band has an area A max = A$e~ 



-At„ 
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determined by the ergodic time r 01g . This is the time 
required for a trajectory to explore the whole accessible 
phase space. The values of r crg and A depend on the 
degree of collimation of the beam of trajectories injected 
into the system*^ For our model, without collimation, 
one has r erg of order unity (one stroboscopic period) and 
Ao ~ (N/M) 2 . The typical transmission band has an 
area Aoe~ XTD which is exponentially smaller than A max 
(since t d = M/2N > r crg ). 

As the position of the lead is moved around, trans- 
mission bands move into and out of the lead. The re- 
sulting fluctuations in the transmission probability Pi ^2 
are dominated by the largest band. Since there is an 
exponentially large number e XTD of typical bands, their 
fluctuations average out. The total area in phase space 
of the lead is A\ ea d = N/M, so we estimate the mean 
squared fluctuations in Pi^2 at 

varfUa ~ (A max /A lcad ) 2 = c{N / M) 2 e~ 2XT ^ , (4.3) 

with c and r crg of order unity. We have tested this func- 
tional dependence numerically for the map l|2.7(l . and 
find a reasonable agreement (see Fig.0J. Both the expo- 
nential dependence on A and the quadratic dependence 
on td = M/2N are consistent with the data. We find 
T org = 0.68 of order unity, as expected. 
Eqs. (|4"2|) and imply 

var G = § + c(iV 4 /M 2 )e~ 2Ar "s . (4.4) 

In Fig. [S] we plot the same data as in Fig. El but now as a 
function of (TV 4 /M 2 )e~ 2Ar ° rs . We see that the functional 
dependence l|4.4|l is approached for large dwell times. 
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FIG. 5: Same data as in Fig. EJrescaled to show the approach 
to a single limiting curve in the large dwell time limit. The 
solid line is calculated from Eq. 14.41 . with the same param- 
eters c = 1.6, T crg = 0.68 as in Fig. 0] 

The quantum fluctuations of RMT dominate over the 
classical fluctuations if iV 2 varPi^2 -C 1. Using the esti- 
mate (I4.3f) . this amounts to the condition 

T crg > max [0, A" 1 ]n(N 2 /M)] = t e (4.5) 



that the ergodic time exceeds the Ehrenfest time. Notice 
that condition (|4.5I) is always satisfied if N 2 < M = 
l/h e g. This agrees with the findings of Ref. 0, that the 
breakdown of RMT starts when N > VM. 



V. CONCLUSIONS 

In summary, we have presented both a fully quan- 
tum mechanical and a semiclassical calculation of the 
quantum-to-classical crossover from universal to non- 
universal conductance fluctuations. The two calculations 
are in very good agreement, without any adjustable pa- 
rameter (compare data points with curves in Fig. |2J. 
We have also given an analytical approximation to the 
numerical data, which allows us to determine the para- 
metric dependence of the crossover. 

We have found that universality of the conductance 
fluctuations requires the ergodic time r org to be larger 
than the Ehrenfest time te- This condition is much 
more stringent than the condition that the dwell time 
td should be larger than te, found previously for uni- 
versality of the shot noise in a quantum dotAiSiii The 
universality of the excitation gap in a quantum dot con- 
nected to a superconductor is also governed by the ratio 
td/te rather than r erg /r£^i2iS, as is the universality 
of the weak localization effectii 2 ^ These two properties 
have in common that they represent ensemble averages, 
rather than sample-to-sample fluctuations. 

We propose that what we have found here for the con- 
ductance is generic for other transport properties: That 
the breakdown of RMT with increasing t e occurs when 
te > td for ensemble averages and when te > r erg for 
the fluctuations. This has immediate experimental conse- 
quences, because it is much easier to violate the condition 
te > Tcrg than the condition te > td- 

To test this proposal, an obvious next step would be to 
determine the ratio of time scales that govern the break- 
down of universality of the fluctuations in the supercon- 
ducting excitation gap. The numerical data in Refs. ^] 
and |2l| was interpreted in terms of the ratio te /td , but 
an alternative description in terms of the ratio te / T eTg 
was not considered. 

One final remark about the distinction between clas- 
sical and quantum fluctuations, explained in Sec. IIIII It 
is possible to suppress the classical fluctuations entirely, 
by varying only the quasi-energy at fixed lead positions. 
In that case we would expect the breakdown of univer- 
sality to be governed by td/te instead of T eTg /TE- Our 
numerical data (Fig. [2J does not show any systematic de- 
viation from RMT, probably because we could not reach 
sufficiently large systems in our simulation. 

Note added: Our final remark above has been criticized 
by Jacquod and Sukhorukov |23. They argue that the 
numerical data of Fig. [21 (and similar data of their own) 
does not show any systematic deviation from RMT be- 
cause quantum fluctuations remain universal if te > td- 
Their argument relies on the assumption that the effec- 
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tive RMT of Ref. holds not only for the classical fluc- 
tuations (as we assume here), but also for the quantum 
fluctuations. The effective RMT says that quantum fluc- 
tuations are due to a number N e g « Ne~ TE / TD of trans- 
mission channels with an RMT distribution. Universal- 
ity of the quantum fluctuations is then guaranteed even 
if iVeff <§; N, as long as N c g is still large compared to 
unity. 

This line of reasoning, if pursued further, contra- 
dicts the established theor} r 12 i 13 of the t e dependence 
of weak localization. RMT says that the weak localiza- 
tion correction 8G = — i is independent of the number of 
channels 3 - 4 . Validity of the effective RMT at the quan- 
tum level would therefore imply that weak localization 



remains universal if te > to, as long as Ne TE / TD ^> 1. 
This contradicts the result SG = L e - TE/TD of Refs. llH 
andd 
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